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Optical resonators provide a powerful tool for testing aspects of Lorentz invariance. Here, we 
present a reanalysis of an experiment where a path asymmetry was created in an optical ring res¬ 
onator by introducing a dielectric prism in one arm. The frequency difference of the two fundamental 
counter-propagating modes was then recorded as the apparatus was orientation-modulated in the 
laboratory. By assuming that the minimal Standard-Model Extension coefficients vanish we are able 
to place bounds on higher-order parity-odd Lorentz-violating coefficients of the Standard-Model Ex¬ 
tension. The results presented in this work set the first constraints on two previously unbounded 
linear combinations of d — 8 parity-odd nonbirefringent nondispersive coefficients of the photon 
sector. 


I. INTRODUCTION 

Lorentz invariance is a fundamental component of the 
Standard Model and General Relativity. Despite the suc¬ 
cesses of both theories they remain incompatible; it is 
generally assumed they are both low-energy approxima¬ 
tions of a single theory that is consistent at the Planck 
scale [If. Various efforts towards identifying a unified 
theory can allow or require Lorentz invariance to be bro¬ 
ken mm- Testing Lorentz invariance thus provides one 
of the few experimental portals for assessing and guiding 
theories of quantum gravity and other unification propo¬ 
sitions. Precision laboratory measurements of Planck- 
scale suppressed effects offer excellent prospects in the 
search for Lorentz-Invariance Violations (LIV) |T. 

The Standard-Model Extension (SME) provides a com¬ 
prehensive framework for analysing, quantifying and 
comparing different experimental tests of LIV [3H6]. The 
SME describes all possible Lorentz and CPT violations 
associated with known particles and fields. Efforts have 
usually focused on the minimal SME sectors, which only 
contain energy-independent operators of renormalizable 
dimension in flat spacetime. In recent years the SME 
has been expanded to include higher-order nonrenormal- 
izable operators, which presents new opportunities for 
experimentation and analysis mm- Bounds have previ¬ 
ously been placed in the photon sector mm, with results 
in other sectors now starting to emerge Dung. Prac¬ 
tically there are situations where one might expect LIV 
to manifest at higher-order, e.g. in some theories with 
noncommutative spacetime coordinates LIV only occurs 
for nonrenormalizable operators [13 E]. 

In the photon sector of the SME operators can be classi¬ 
fied into different groups that describe their effects on 
standard electrodynamics. Astrophysical observations 
have tightly constrained birefringent and vacuum disper¬ 
sive effects with sensitivities far beyond the reach of ter- 
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restrial tests. This provides physical motivation to the 
study of the camouflage coefficients, which to leading or¬ 
der are nonbirefringent and nondispersive. The camou¬ 
flage coefficients are CPT invariant, have even dimension 
d > 4, and are best constrained via precision electromag¬ 
netic resonant-cavity tests, such as the modern descen¬ 
dants of the Michelson-Morley experiment. 

Constraints have previously been set on combinations 
of d = 6 and d = 8 camouflage coefficients using data 
from a parity-even microwave cavity experiment [9| and 
a parity-odd optical cavity experiment m- As the sen¬ 
sitivity to the camouflage coefficients scales with the fre¬ 
quency of the photon (roughly as z/ d-4 ), optical cavities 
are very well suited for testing these effects [3- For the 
parity-odd camouflage coefficients there are three d = 6 
coefficients and 13 d = 8 coefficients, with three lin¬ 
ear combinations currently unbounded. Here we analyse 
an orientation modulated parity-odd asymmetric opti¬ 
cal ring resonator experiment to set bounds on higher- 
order parity-odd camouflage coefficients, including the 
first constraints on two of the three remaining unbounded 
d = 8 combinations. 


II. EXPERIMENTAL SETUP 

A detailed overview of the experiment is provided 
in [T6| and [15] . As seen in Fig. [T] a dielectric prism with 
index of refraction n = 1.44 is placed in one arm of an 
optical ring resonator to create a path asymmetry. A 
1064 nm laser is split and each path is independently 
locked to a fundamental counter-propagating mode. The 
first laser path is frequency shifted by an Acousto-Optic 
Modulator (AOM), introducing a 160 MHz offset, and 
then mode-locked by applying correction signals to 
the laser. The second path is locked to the counter- 
propagating mode via an AOM; the correction signal 
contains the frequency difference (typically ^100 Hz) of 
the two modes imprinted on the 160 MHz offset. This 
was recorded directly with a frequency counter. 
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FIG. 1. Schematic of the parity-odd asymmetric optical ring resonator experiment, reproduced from m- This only shows the 
control systems and setup for one of the counter-propagating modes, there is a nominally identical setup for the second mode. 


The apparatus sits upon a rotation platform. The beat 
frequency was recorded for 10 minutes in a stationary 
position, then the experiment was rotated by +180° 
and the beat frequency recorded for 10 minutes while 
stationary. The experiment was then rotated back by 
-180° to the initial position and the cycle repeated. 
The result of this is that the final dataset is in essence 
stationary, but with leading-order noise processes shifted 
to a different frequency. Previously, measurements were 
made for 50 days and used to place bounds in the mini¬ 
mal SME [T5j . We use the same dataset for the following 
analysis, but assume that the minimal SME terms vanish. 


III. EXPERIMENTAL SENSITIVITY TO 
HIGHER-ORDER CAMOUFLAGE 
COEFFICIENTS 


A complete derivation of the camouflage coefficients 
is provided in [|7]. In short, the coefficients associated 
with CPT-even differential operators from the arbitrary 
Lorentz and CPT violating Lagrange density for photon 
propagation are decomposed via spin-weighted spherical 
harmonics. This then allows for groups of coefficients 
associated with different physical properties to be col¬ 
lected together. The coefficients from operators with no 
leading-order birefringence are denoted where 

d is the mass dimension, 0 is the spin weight, E = (—1 ) J is 
the parity, n is the wavelength dependence, j is the total 
angular momentum and m is the ^-component of angu¬ 
lar momentum. The camouflage coefficients are then the 
subset of coefficients that are not associated 

with leading-order vacuum dispersion, 


, (d)x(0S) _ /-.(d)x(O E) 
\ U F Jnjm V C F Jnjm 


(7 ,(d)A0E) 
v F J(n—2)jm‘ 


Here, we restrict our attention to d=6,8, with 0 < n <(d- 
4), j = n, n — 2,..., >0 and \m\ < j. For coefficients with 
m/0 there are real and imaginary components. 


The measured fractional beat frequency between the two 
counter-propagating modes of Fig. [l] is 

^ beat — A rnrn /e im ^+ irn ' UJ ® T ® , (1) 

mm' 

where </> is the angle between the laboratory-frame x axis 
and geographical south, c is the Earth’s sidereal rota¬ 
tion rate and T 0 is the sidereal time since the last align¬ 
ment of the experiment with the Sun-Centred Celestial 
Equatorial Frame (SCCEF), the conventional choice of 
reference frame for analysis of such experiments within 
the SME [5]. The A mm / factors contain linear combi¬ 
nations of Lorentz-violating coefficients that a given ex¬ 
periment is sensitive to. They satisfy the relationship 
A^m' = an d can be calculated from 

A mm , = Y A(2) 

dnj 

Here AM^ a6 is an experiment-dependent constant 
that considers the difference between the two counter- 
propagating modes of Fig. HI d are the little Wigner 
matrices and % is the co-latitude of the experiment. A 
detailed guide describing these factors and how to calcu¬ 
late them is provided in section IV of [8]. 

As is standard for the analysis of turntable experi¬ 
ments 0E! we can model the data as 

^ beat = ^ Cm cos (me/)) + S m sin (ra</>), 

m> 0 

( 3 ) 

Cm = Y C ™ m ' C0S ( m ' W © T ©) + Cmm' sin (m'w® 1 ©). 

m'> 0 

( 4 ) 

s m= Y S ™rn’ cos (™'^© T ©) + S^ m , sin (m'w®?©), 

m'> 0 

( 5 ) 
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TABLE I. Sensitivities of sidereal amplitudes to d = 6 
(10 -18 GeV 2 ) and d = 8 (10 -36 GeV 4 ) camouflage coeffi¬ 
cients. Amplitudes are denoted as (7®, for cosine and S'®, for 
sine; m is an integer harmonic of c j®T® as discussed in the 
main text. 


Amplitude 

Sensitivity 

ZD 

II 

73 


cf 

-0.37 Im(40ui 

sf 

-0.37 Re^lOm 

d = 8 


Cf 

Im(c?!’)iii - 10.76 Im(7?^) 3 n + 6.00 Im(7?^) 33 i 

Sf 

Re(40in - 10.76 Re(c l |') 3 n + 6.00 Re(40 33 i 

cf 

3.25 Im(^) 33 2 

sf 

3.25 Re (?|) 3 32 

cf 

10.44 Im(c’|') 333 

Sf 

10.44 Re(40 33 3 


with the cosine and sine amplitudes of Eqs. 0 and 0 
given by 

^mm' ^VmVm' Re( A mm / T 

Cmm' = — — A m (_ m /)), 

^rnm’ Im(A mm / T A m (_ m q), (6) 

^mm' = — 2Re(A mm / — A m (_ m q), (7) 

with 770 = 1/2 and 7? m = 1 for all other values. Despite 
the experiment being mounted on a turntable the dataset 
is effectively stationary with the experiment x-axis per¬ 
manently aligned parallel to East-West, thereby setting 
4> = 7t/ 2. This sets the cosine component of Eq. 0 to 
zero, so experimental access to Lorentz violating cam¬ 
ouflage coefficients is restricted to Eqs. © and 0 . In 
practical terms the beat frequency will be fit to a model 
containing a linear offset and quadrature components of 
sidereal harmonics (m' = 1,2,3), 

= a + ^2 C m‘ C0S + S m' sin (V^©r©)- 

m'> 0 

(8) 

Using equations 0-0 we can explicitly determine the 
sensitivity of the experiment to higher-order coefficients; 
results are summarized in Table [l| The amplitudes Cf 
and Sf provide access to 2 combinations of d = 8 coeffi¬ 
cients that have yet to be bounded. 


IV. RESULTS 

Data analysis proceeds in a similar manner to previous 
work mum. The frequency difference between the two 
counter-propagating modes is averaged while the appara¬ 
tus is stationary (~10 minute period). This collection of 
averaged data points provides us with the Vbeat term in 


TABLE II. Measured amplitudes of sine and cosine compo¬ 
nents of 7 ?t/u;®T®, where m is denoted in the amplitude co¬ 
efficient subscript. All values are given x 10 -14 . Errors are 
2(7. 


Amplitude 

Measurement 

cf 

5.60 

± 

3.96 

sf 

5.73 

± 

3.98 

cf 

-1.17 

± 

3.96 

sf 

6.56 

± 

3.98 

cf 

-0.44 

± 

3.98 

sf 

-2.90 

± 

3.96 


Eq. |3]). Ordinary Least Squares regression is used to si¬ 
multaneously fit to an offset and quadrature components 
of modulation of the form a;®, 2 cj® and 3 cj®. The fit is 
performed over subsets of the beat frequency, each com¬ 
prising 120 points, which is equivalent to one day. This 
minimizes the standard error of the fit. The magnitude of 
the error for each frequency of modulation is computed, 

y^rrL© + 5 an d an y day °f measurement where the 

magnitude of the error is more than 3cr from the mean 
is discarded. This resulted in the removal of 1 day from 
the dataset, leaving 50 values for each amplitude which 
are then combined via an error-weighted average, 



y A 
^ E 2 

y 4 

^ E 2 

(9) 

1 

V 1 ’ 

^ E 2 

(10) 


where A is the amplitude and E is the corresponding 
standard error. Final values with 2cr errors are presented 

in Table HI 

Following the convention established in m we also take 
the error-weighted average of all cj® and 2 cj® modula¬ 
tion amplitudes to determine an overall constraint on a 
shift in the fractional frequency, A v/v <4.2±2.0xl0 -14 
(95% C.I.). The purpose of such a value is to provide a 
single number to represent the fundamental sensitivity of 
the experiment. 

The Sf amplitude is statistically significant at a 3.do- 
level, however this significance is not constant over the 
50 day run of the experiment and only manifests during 
a ~10 day subset. By plotting the quadrature amplitude 
fits for each modulation frequency against each other we 
do not see any persistent signal congruent with a viola¬ 
tion of Lorentz invariance (Fig s. |2| [ 3 ] and |4| . 

Now we use Table [T] and Table [Tlfto place constraints on 
higher-order Lorentz-violating coefficients of the SME; 
results are displayed in Table EH Our constraints for 
d = 6 and d = 8 coefficients already bounded by previous 
work m are 2 orders of magnitude larger, which is con¬ 
sistent with the differences in absolute sensitivity, dura¬ 
tion of measurements and rotation rate between the two 
experiments. Using our first bounds for two of the linear 
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combinations of d = 8 coefficients we are able to infer the 
first constraints for Re( 0 ^) 311 , Im(c , ^) 3 n, Re(c , ^)m 
and Im(c 1 ^)in. Unfortunately, due to the East-West 
alignment of the apparatus, the experiment cannot be 
used to place a limit on the one remaining d = 8 lin¬ 
ear combination, ( 0^)310 - 0.093 However, all 

of the other parity-odd nonrenormalizable higher-order 
d = 6 and d = 8 coefficients of the photon sector of the 
SME now have at least one bound constraining them. 
Recent work from a parity-even microwave cavity test of 
Lorentz invariance achieved a fractional frequency sensi¬ 
tivity of 10 -18 [17]. A planned co-rotating optical and 
microwave cavity setup with an estimated sensitivity be¬ 
yond this pa uni will enable significant improvements 


to both minimal and parity-even and parity-odd higher- 
order nonrenormalizable SME coefficient bounds. 
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TABLE III. Bounds on higher-order SME camouflage coefficients, d = 6 coefficients have units GeV 2 and d — 8 coefficients 
have units GeV -4 . Comparisons are given to previous work m , combining results from both experiments allows us to infer 
some individual bounds. Errors are 2 a. 


Coefficient 

This Work 

Previous Work ^10 

d = 6 



Re(c , |)in 

(-1.55 ± 1.08) x 10® 

(-0.8 ± 2.2) x 10 3 

Im(40in 

(-1.51 ± 1.08) x 10® 

(-0.6 ± 2.2) x 10 3 

00 

II 

"0 



(c?f)3io ■ 0.020 (c'l’)iio 
Re(E^)3ii - 0.020 Re(c l |’)in 

Im(c , |’) 3 n - 0.020 Im(c , | 1 )in 

Re(^) 3 u - 0.093 Re(^) m 

(-5.33 ± 3.70) x 10 21 

(-0.2 ± 3.8) x 10 19 
(1.4 ± 2.6) x 10 19 
(0.1 ± 2.6) x 10 19 

Im(c l |’) 3 n - 0.093 Im(E|)m 

(-5.21 ± 3.68) x 10 21 

- 

Re(c l |^) 33 2 

(20.19 ± 12.24) x 10 21 

(2.2 ± 2.6) x 10 19 

Im(c l |^) 33 2 

(-3.60 ± 12.24) x 10 21 

(0.2 ± 2.6) x 10 19 

Re(c'|r)333 

(-2.78 ± 3.80) x 10 21 

(-0.1 ± 3.2) x 10 19 

Im(E^) 33 3 

(-4.22 ± 3.82) x 10 21 

(-0.1 ± 3.2) x 10 19 

Inferred 



Re^fOsn 

(1.48 ± 0.98) x 10 21 

- 

Im(c , ^) 3 n 

(1.43 ± 0.98) x 10 21 

- 

Re(^)m 

(7.32 ± 5.04) x 10 22 

- 


(7.14 ± 5.00) x 10 22 

- 
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FIG. 2. Fitted amplitude of fractional frequency shift for 
the cu® modulation with la error bars. Each data point is 
obtained by fitting to ~1 day of data, as described in the 
main text. 



-2E-12 -IE-12 0 IE-12 2E-12 


Cosine Amplitude 

FIG. 4. Fitted amplitude of fractional frequency shift for 
the 3 cj® modulation with la error bars. Each data point is 
obtained by fitting to ~1 day of data, as described in the main 
text. 



FIG. 3. Fitted amplitude of fractional frequency shift for 
the 2 cj® modulation with la error bars. Each data point is 
obtained by fitting to day of data, as described in the main 
text. 








































































































